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We use the random phase approximation to systematically describe the total muon capture 
, rates on all nuclei where they have been measured. We reproduce the experimental values 

' on these nuclei to better than 15% accuracy using the free nucleon weak form factors and 

residual interactions with a mild A dependency. The isospin dependence and the effects 
associated with shell closures are fairly well reproduced as well. However, the calculated 
rates for the same residual interactions would be significantly lower than the data if the 
in-medium quenching of the axial-vector coupling constant is employed to other than the 
VO ■ true Gamow- Teller amplitudes. Our calculation thus suggests that no quenching is needed 

in the description of semileptonic weak processes involving higher multipole transitions and 
momentum transfer ~ m M , with obvious importance to analogous weak processes. 
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PACS numbers: 24.30. Cz, 23.40.-s, 23.40. He 



INTRODUCTION 



The capture of a negative muon from the atomic Is orbit, 

5-H ' 



l r + (z ) N)^v li + (z-i,N + iy (i) 



is a semileptonic weak process that has been studied for a longtime (see, e.g., the recent review |l| or the 
earlier one by Walecka 

2] 

and the classic by Mukhopadhyay 3] and the earlier references therein). The 
total capture rate has been measured for many stable nuclei; in some cases the capture rates on separated 
isotopes have been also determined 0|. 

The nuclear response in muon capture is governed by the momentum transfer of the order of the muon 
mass. The phase space and the nuclear response favor lower nuclear excitation energies, thus the nuclear 
states in the giant resonance region dominate. Since the experimental data are quite accurate, and the 
theoretical techniques of evaluating the nuclear response in the relevant regime are well developed, it is 
worthwhile to see to what extent the capture rates are understood globally. Such a comparison may be 
viewed as a general test of our ability to describe semileptonic weak charged-current reactions with q ~ 
over a large range of nuclei, where q is the momentum transfer and m M = 105.6 MeV is the muon mass. 

The present work represents a first fully comprehensive theoretical evaluation of the total muon capture 
rate over the full range of nuclei where the experimental data are available. Previously, the muon capture 
rates for selected nuclei encompassing, however, a broad range of atomic charges, were calculated in Ref. 
[jj. That paper was devoted mostly to the description of the radiative muon capture, and the total muon 
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capture rates were a byproduct with only a limited agreement with the data. More along the lines of the 
present approach, Ref. p and 0] used a Hatree-Fock random phase approximation method and obtained 
good agreement with experiment, however only for a limited selection of nuclei. The local Fermi gas model 
was used successfully for the evaluation of the muon capture rate in selected nuclei in Ref. 0], and more 
recently in Ref. 

The present work is an extension of previous papers devoted to this issue 

H3, EH In Ref. the 

capture rates for 12 C, 16 and 40 Ca were evaluated using the continuum random phase approximation, 
and a very good agreement with the total rate was obtained. However, the residual interaction employed 
in ^(j] was adjusted to describe other observables in the cases of 12 C and 16 0. Moreover, it was necessary 
to quench the Gamow- Teller like (GT) partial capture rate leading to the 1 + ground state of 12 B. In the 
later Ref. heavier nuclei with N > Z, 44 ' 48 Ca, 56 Fe, 90 Zr, and 208 Pb were also included with a similar 
success. In Ref. [Q.] it has also been shown that for the calculation of muon capture rates the standard 
random phase approximation (SRPA) is essentially equivalent to the more computationally demanding 
continuum random phase approximation. Thus, the SRPA method is also used in the present study. In 
Ref. 12] the SRPA approach has been used to study the muon capture rates to a long chain of calcium 
and tin isotopes. 

One of the important issues when evaluating the response of nuclei to weak probes of relatively low 
energies is the problem of quenching of the corresponding strength. The evidence for quenching comes 
primarily from the analysis of the beta decay of the (sd) shell [l^l as well as (p, f) shell nuclei. In 
addition, the interpretation of the forward angle (p, n) and (n,p) charge-exchange reactions jlil. Hil 
leads to the same conclusion. All such evidence, so far, is restricted to the GT strength and relatively low 
excitation energies. A convenient and customary way to account for this quenching is to use an 'effective' 
axial- vector coupling constant gA, reducing it from its nominal value of gA = 1-26 to gA ~ 1. 

The evaluation of the muon capture rate, reported here, suggests that the quenching of gA is not needed 
to describe these data. (That conclusion was already reached in Refs.EJ As stressed above, the mean 
excitation energy in muon capture is in the region of giant resonances of about 15 MeV (slowly decreasing 
with A or Z), and the GT-like operators contribute very little in heavier nuclei where the neutrons and 
protons are in different oscillator shells. In lighter nuclei, for N and Z less than 40, the GT strength 
contributes, and is concentrated at low energy. Thus, in agreement with the evidence mentioned above, 
we quench this, and only this, part of the transition strength by a common factor (0.8) 2 = 0.64 [l3| . 

The present work, and the evaluation of the muon capture in general, makes it possible to extend the 
study of quenching to higher multipoles, and correspondingly to higher nuclear excitation energies. Such 
processes, typically, depend primarily on the positions of the corresponding giant resonances and on the 
overall strength. Our conclusions, therefore, show that the SRPA method involving correlated particle-hole 
excitations, is capable of describing the inclusive semileptonic processes with momentum transfer q of order 
~ quite well. This is an important conclusion, applicable to a variety of practically important subjects, 
e.g. detection of supernova neutrinos or evaluation of the nuclear matrix elements for neutrinoless double 
beta decay. 

The challenge of evaluating the muon capture rate in a wide variety of nuclei made it necessary to include 
several effects that were not usually included in analogous calculations. Since we needed to describe the 
bound muon in the Is orbit well, we went beyond the usual calculation of the muon density at the site 
of the nucleus. First, we solved the Dirac equation in the field of the finite size nucleus numerically. We 
then used its wave function, taking into account that it is not constant between the origin and the nuclear 
surface. For high Z values the muon is relativistic, and the 'small' P\i2 component of its wave function 
is nonnegligible. As explained below, we used here (for the first time) the additional transition matrix 
elements associated with that component. 

Since our goal is to describe muon capture in all nuclei (except the very light ones) we have to desribe, 
at least crudely, effects associated with the partial filling of the single-particle subshells for nuclei that do 
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not have magic numbers of protons and/or neutrons. As described in the next section, we describe these 
effect by taking into account the smearing of the proton and neutron Fermi levels caused by pairing and 
deformation. It appears that this simplified treatment of complicated correlations, including those caused 
by deformation, is sufficient for our purpose. 

II. METHOD AND PARAMETERS 

In this calculation we have used the standard RPA model to describe the nuclear excitations. In a previous 
work ([lljp this model was shown to be just as good as the computationally more involved continuum 
RPA. As residual interaction we use the phenomenological Landau-Migdal force. For low mass nuclei the 
parameters for the force were taken from Ref. This choice was shown to be accurate in For 

muon capture the most important term in the Landau-Migdal force is the spin-isospin coupling constant 
g'. In [18( the value g'=0.7 is recommended, however, in a recent review [l9| g'=0.96 is used for heavy 
nuclei. To accommodate this variation, we use an interpolation formula with a mild A dependency: 

g' = Cl + c 2 A 1 / 3 (2) 

where the constants c\ and c 2 are fitted to yield g'=0.7 in 16 and g'=0.96 in 208 Pb. We note that the 
change in the total capture rate in going from g'=0.7 to g'=0.96 is less than 10%. 

To get a basis of single-particle states we diagonalize a Woods-Saxon potential (WSP) in a harmonic 
oscillator basis of more than 8 major shells, thus enabling us to always have an excess of 2fuv of valence space 
above the Fermi level for both protons and neutrons. As parameters of the WSP we use 

#0 = 1-2 * A 1 / 3 

fm for the radius and a = 0.53 fm for the diffuseness. The spin-orbit term is given as the derivative of the 
WSP times a strength V so . Here we have simply used a fixed value of V so = —8.95 through-out, initially 
checking that other choices did not significantly effect the total capture rate. To find the overall strength 
of the WSP, we fixed the last proton and neutron particle energies to experimentally known masses. More 
specifically, for a nucleus (A,Z) we found the energy of the last proton level from the proton separation 
energy S p in (A,Z). For the last neutron level we used the neutron separation energy S n , but this time in 
the daughter (A,Z-1). 

In order to be able to handle open-shell nuclei we have previously used a simple scheme where partial 
occupancies was treated by multiplying the open level matrix elements by occupation numbers correspond- 
ing to an independent particle model |2jJ. In this work we have attempted to improve on this treatment 
by solving the standard BCS equations to determine the occupation numbers. Following Ref.ji3] the 
occupation probabilities are given by 

where e/% are the single-particle energies and the chemical potential [i is fixed by the condition N = 
Y^k v k(2jk + !)• The pairing gap A is obtained by the procedure descibed in Ref. [23j. 

The formalism used to evaluate the total muon capture rate is that of Ref. As mentioned in the 
Introduction, we treat the muon wavefunction by solving the Dirac equation in the extended charge of the 
nucleus, which is assumed to be of the Woods-Saxon form with the same paramteres as above. For nuclei 
with large values of Z, the atomic binding energy becomes a significant fraction of the muon rest mass 
and the small component of the Dirac bi-spinor may not be negligble in this range. We therefore explicitly 
include all terms containing both large and small components in our transition operators. An outline of 
the complications arising from this is given in the Appendix. 
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III. RESULTS 
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FIG. 1: (color online) Comparison between the measured total muon capture rates denoted by squares □, the 
calculated rates with all corrections (empty circles O), an d the calculated rates without the BCS and relativistic 
corrections (diamonds ). The insert, in larger scale, shows the same results for light nuclei. When the measurements 
are for the natural abundance of a given element, the calculation represent the corresponding combination of the 
individual isotopes. 



Our calculated total muon capture rates for all nuclei for which measured values exist are shown in Figure^ 
One sees that the overall agreement is quite good. With the exceptions at Z = 74 and 80 the calculations 
reproduce the experimental values to 15% or better. In Figure El we provide a ratio plot where the degree 
of agreement is better seen. In addition, in Table [I] we collect all our calculated capture rates, including 
the results for individual isotopes. 

In Fig. ^the total capture rates are also compared with values obtained from calculations where the 
BCS occupancy and the relativistic corrections are turned off. From this comparison it is clear that these 
are only small corrections, providing justification for the calculation done in Ref. 

0,1 

Til . Eol ] . At low Z one 

especially notices the good reproduction of the distinct dips in the rates above the magic numbers Z = 20 
and 28. For Z = 50 the same trends is visible in both calculation and experiment, but the calculations 
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overshoot the experimental values somewhat. Just above Z = 60, where the N = 82 closed neutron shell 
comes into play, we also overestimate the capture rates. This continues into the region with Z = 74 and 
80, and also, to a lesser extent, to the doubly-magic nucleus 208 Pb. One should remember that some of the 
nuclei above Z = 50 and below Z = 82 are deformed and thus have different single-particle structures than 
the ones given by our spherical mean-field model, thus a perfect agreement should not be expected. The 
fact that the calculated values are again approaching experiment at Th,U and Pu is likely a consequence 
of the same dip after the magic shell closure that was seen at lower Z also (these trends were also noted 
in Ref. 2J|, see in particular figure 3 of that reference). 

As stated above, we use the unquenched value for the axial-vector coupling constant for all multipole 
operators, except for the true Gamow- Teller transition. For most of the light- and medium-mass nuclei, 
A = 1 — (dipole-like) and A = 2 + (quadrupole-like) transitions dominate. However, for 208 Pb and the 
heavier nuclei, A = 1 + transitions contribute significantly to the total capture rate. For these nuclei, the 
neutron excess is already so large that, in the simple independent particle model, 2 major oscillator shells 
have to be overcome when changing a proton into a neutron in the muon capture process. Thus, for these 
excitations the A = 1 + multipole transition corresponds to a 2hw mode and not to a (0 fvuS) Gamow- Teller 
transition; the contribution of the latter to the rate in the heavy nuclei vanishes in our calculations. We 
have not renormalized the axial-vector coupling constant for such 2 huj 1 + transitions, supported by the 
good agreement of the the calculated capture rates with the measured results for the heavy nuclei. 

An important test of the ability of a nuclear structure model to reproduce the data is the dependence 
of the muon capture rate on the number of neutrons for a fixed nuclear charge Z. Typically, the rate 
decreases with increasing N (or A) as subsequently more neutron levels are getting blocked. This effect is 
incorporated into the well-known Primakoff parametrization by its (N-Z) dependence [2^]. As examples, 
table [I] include three isotope chains of nuclei (Ca, Cr, and Ni), where total capture rates for individual 
isotopes have been measured. One sees that the isotope dependence is well reproduced by our calculations. 
Analogous calculations were performed in Ref. 0,0. 



IV. CONCLUSION 

The present analysis shows that the standard random phase approximation method is capable of describing 
quite well the total pT capture rates for essentially all stable nuclei. The dependence of the capture rate 
on the isospin, or neutron excess, the so-called Primakoff rule [2||, is also fairly well, albeit not perfectly, 
reproduced. Our calculation even describes the rather subtle effects of shell closures when considering the 
dependence of the capture rate on Z and/or A. There is no indication of the need to apply any quenching 
to the operators responsible for the muon capture, in particular those involving single-particle transitions 
from one oscillator shell to another, i.e. other than those involving 0/zu; spin and isospin changing operators. 

Given the task of describing the capture in a variety of nuclei, including those with high charge Z 
and nuclei with unfilled shells, it became necessary to consider several effects that have not been typically 
included previously. One of them is a fully relativistic treatment of the muon bound state, including the 
effects associated with the 'small' P\/2 component of its wave function. Another one is the effect of the 
smearing of the Fermi level (both for protons and neutrons) in nuclei that have nonmagic Z or N numbers. 
Even though the corresponding corrections are not very large, they contribute noticeably to the overall 
good agreement between the experimental data and our calculated values. 

Our findings then can be used as guidance in the evaluation of a wide variety of semileptonic weak 
processes on nuclei with similar momentum transfer. 
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FIG. 2: (color online) Ratios of the calculated and measured total muon capture rates vs. the atomic number. 
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VI. APPENDIX 



In this appendix we give an outline of terms arising from the inclusion of the fully relativistic treatment 
of the bound muon wave function, i.e., of both the large and small components. The muon wavefunction 
we use has the form 



IJ jm 



ig(r)ttji m (r) 
-f{r)VL jVm (f) 



(4) 



where the radial functions satisfy the equations 




_l±2 E + m^-V{r 

-E + m^ + V{r) 




(5) 



where g(r),f(r) are the large and small components, respectively. Here 

(6) 



-{1 + 1) ,3 = 1 + 1 
I , 3 = l~\ 



and 

ttj,l,m(r) = ^ (l,l/2J\mi,m a ,m)Y ljmi (f)xm a (7) 



mi,m s 



These equations are entirely general and can be found in e.g. Ref. |2£(]. Since we assume that the muon 
is captured from the atomic Is orbit we have I = and V = 1. Since V is non-zero, we can no longer just 
multiply the wavefunction with the irreducible nuclear operators and obtain good total angular momentum. 
If the positive z-axis is chosen along the direction of the out-going neutrino then its wavefunction becomes 

( X- \ 
ip = N I I exp(ip • x) 

Here X- are the usual spin down Pauli two-spinors and N is a normalization given in [2^]. This wavefunction 
makes the neutrino purely left-handed as the standard model prescribes. 

The approach used in 0] neglects the small component in the muon wave function and expands the 
neutrino plane wave in multipoles. As the large component has I = 0, angular momentum coupling of the 
muon wave function and the multipole operators is quite straightforward. This is, however, no longer the 
case if the small component with orbital angular momentum I = 1 is considered, implying the need for a 
cumbersome recoupling of angular momenta to regain tensor operators that can be applied in the nuclear 
Hilbert space. This results in a more complicated expression for the weak Hamiltonian governing muon 
capture with several new terms. In the notation used in |2| the Hamiltonian with all terms from both 



components can be written as 



oo 

[£V&[J](-i) J {i6 m> - 1/2 {M'j fi 



n = 

_2G F cos9 c N* 

~ ^ ,=o 
+a (j-l,J,m + \) Ti (j-l,J,m+ \) 

+q ( J + 1, J, m + ±) Ti ( J + 1, J, m + \) 
-i(3 + (J, J, m)T 2 (J, J,m + \) 
-i(3+(J + 1, J, m)T 2 (J + 1, J, m + ±) 
—i(3 + (J + 1, J + 2, m)T 2 ( J + 1, J + 2, to + ±) 
-i(3-(J — 1, J — 2,m)%(J — 1, J — 2,m + |) 
-i/3-(J- 1, J,m)T 3 (J- 1, J,m + ±) 
-i/3-(J, J,m)T 3 (J, J,m + i)} 

oo 

+ £ V4^[J](-i) J {i<W/2 { J# - jJt} 
j=i 

— 5(J - 1, J — 1, to)T 4 ( J — 1, J — l,m + |) 

— <5(J — 1, J — 1, m)T 4 ( J — 1, J — l,m + |) 

—5(J, J, m)Ti(J, J, to + ^) 

-<5(J, J + 1, m)T 4 (J, J + 1, m + §) 

— <5(J + 1, J + 1, to)T 4 ( J + 1, J + 1, to + ±) 

+i7] + (J, J, m)T 2 (J, J,m + \) 

+ii] + {J + 1, J, m)T 2 (J + 1, J,m + |) 

+ii] + (J + 1, J + 2, m)T 2 (J + 1, J + 2, to + ±) 

— irj-(J — 1, J, m)%(J — 1, J, m + |) 

-irj-(J, J,m)T 3 (J, J,m + \) 

-ir,_{J - 1, J - 2, m)T 3 ( J - 1, J - 2, m + 
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Here we have defined the tensor operators in the nuclear Hilbert space as 



M J>M 

J 'mag 
J,M 



n-'el 
J J,M 



d 3 xg(r)Y 0fi jj(Kx)Yj :M Jo 



d 3 £g(r)y ,oV {jj{Kx)Y JM ) Jo 



A 3 xg{r)Y 0fi jj{Kx)yf )Jtl ■ J 



M 



- / d* J %(r)y 0j0 V A [jj{Kx)yfji) ■ J 



p,p) 
%(n,p,j,ii) 

%(i,p,j,p) 



d xf(r)j p (Kx)Y 7!ll Jo 
d^xf(r)jj+i(Kx)y^ pl ■ J 
d 3 xf(r)jj^(Kx)y!; >P;1 -J 
d 3 xf(r)jj(Kx)y^ pA ■ J, 



where Y are the spherical harmonics and y are the vector harmonics. The first four operators are those 
involving the large component and are identical with the ones given in Q|. Their tensor character is (J, M). 
The last four are new operators, i.e. they were ignored in |2(, involving the small component of the muon 
wavefunction. Their tensor character is (7, p). The other indices of the new operators identify terms which 
originate in the multipole expansion to produce the correct spherical Bessel function in the integrals. The 
constants appearing in the Hamiltonian above are given by the rather lengthy expressions 



a(j,p,m 
P+(l,P,J,m 

P-(p(,P,J,m 
5(7, p, J, m 

V+(jiP,J, m 
r/_(7, p, J, m 



3/2+m 
3 



^M( p i 7 |000) <pl 7 |0m + im + i; 



2^^x/jTT[J+l](lJ+lp|000) 
x (U7|m + i0m + !) W(l J + l~fl;pJ) 

^T^^Ji^JV - i](i J - 1 pIooo) 

<U7|m + ±0m + \)W (I J - ljl;pJ) 
X (Uy\m- \lm + \) W(lJ^l;pJ) 

3J ^^Vj[j + i](ij + 1p\ooo) 

x (U7|m- ±lm + \)W (1 J + IjhpJ) 



^^VjTT[j-i]<ij-Mooo) 

x (U-f\m - \ 1 m + |) W (U - 1 7 1; p J) 



Here we have repeatedly used the standard notation [J] = \/2J + 1. All conventions for the Clebsch-Gordan 
coefficients and W-symbols are those of (29] ]. 

In the derivation of the Hamiltonian given above we have used various selection rules for the Clebsch- 
Gordan coefficients and the W-symbols. Note that some of the terms vanish at low angular momenta, since 



10 



7 and p must be positive for any term to contribute. The quantity m in the above expressions corresponds 
to the spin projection of the muon. Since we consider unpolarized muons we must average over the two 
values m = ±1/2. 

For completeness we list here the relevant nuclear currents for the muon capture process. These are 



3=1 
A 



J v (x) 



Ge 
2Mi 



3=1 



2M 



3=1 



2Mi 



3=1 



A 



^V-^T + (i)a(i)^)(^-^) 

3=1 



2M 



J A (x) = G A ^r+0V0')* (3) (^-%) 

3=1 

here M is the nucleon mass, Ge and Gm are the Sachs nucleon form factors and Ga is the axial form 
factor. We note that that the usual Fermi and Gamow- Teller transition operators are recovered in the 
q — > limit as the following multipole components (see, e.g. |3pl |) 



1 A 
M 0fi = -7=GyVr + (j) 

V47T _ 



£i,M - -7=T\,m 



3=1 

1 n-el i 

v 7 ! XM = 



12:; 



G^^T + 0')cr liA4 -(i) 



were Gy is the g — ► limit of the vector coupling form factor, which is often denoted by -Fi(O) in the 
literature. 

To get the final expression for the total rate one must now evaluate the absolute squared matrix element 
of the Hamiltonian in the initial and final nuclear states and multiply it by the two-body phase space factor 
given in It is advantegeous to group together the components with like tensor order to better control 
the interference of operators arising from the large and small components. 
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Nuc 


exp 


calc 


Nuc 


exp 


calc 


Nuc 


exp 


calc 


1 9 S~~i 

c 


0.039 


0.042 


1 fi/~v 
IbQ 


0.103 


0.104 


1 8 /~\ 
18Q 


0.088 


0.089 




0.204 


0.206 


94t\ t 


0.484 


0.454 


28 o- 
9 Si 


0.871 


0.823 


S9 n 


1.352 


1.269 


40 a 
4U Ar 


1.355 


1.345 


4n ft 


2.557 


2.379 


44 Ca 


1.793 


1.946 


48 /~i 

48 Ca 


1.214 1 


1.455 


48m* 

48 Ti 


2.590 


2.214 


n Q t Cr 


3.472 


3.101 


0U Cr 


3.825 


3.451 


^2 n 


3.452 


3.085 


54 Cr 


3.057 


3.024 




4.411 


4.457 


nat Ni 


5.932 


6.004 


08 Ni 


6.110 


6.230 


bu Ni 


5.560 


5.563 


fi2l^T- 

b ^Ni 


4.720 


4.939 


nat Zn 


5.834 


5.235 


b4 Zn 




5.735 


bb Zn 




4.976 


bs Zn 




4.328 


natQ Q 


5.569 


5.317 


7n si 

7U Ge 




5.948 


' 2 Ge 




5.311 


74 /"i 

' 4 Ge 




4.970 


natgg 


5.681 


5.588 


' s Se 




6.023 


su Se 




5.485 


82 n 




5.024 




7.020 


7.529 


sbg r 




8.225 


S8 Sr 


6.610 


7.445 


7 1 0, t ^ j-. 


8.660 


8.897 


yu Zr 




8.974 


Q9 r~i 

y2 Zr 




9.254 


Q4r7 

y4 Zr 




8.317 


""'Mo 


9.614 


10.33 


Q9t\ it 




10.80 


Q4 tv x 

y4 Mo 




11.01 


Mo 




10.04 


y »Mo 




9.153 




10.00 


11.00 


iU4 Pd 




12.71 


1UD Pd 




11.44 


iU8 Pd 




10.44 


liu Pd 




9.607 


™ at Cd 


10.61 


11.46 


linn i 

iiU Cd 




12.58 


1 1 9 s~i i 




11.51 


ii4 Cd 




11.21 


iiD Cd 




10.44 


natg n 


10.44 


11.95 


1 1 fin 




13.08 


1 1 H n 




12.35 


1 2fln 




11.64 


1 22 n 

^Sn 




10.82 


1 24 o 




10.15 




9.270 


9.523 


1 2firTH 

iZb Te 




10.20 


1 2Km 




9.639 


1 30m 




9.043 


not Ba 


9.940 


10.80 


i3b Ba 




11.45 


1 38"T» 




10.73 


" a *Ce 


11.60 


12.44 


i 4n s i 

140 Ce 




12.38 


1 49 r t 

142 Ce 




12.95 


nat Nd 


12.50 


13.70 


1 49 TVT 1 

i4 ^Nd 




13.67 


144 Nd 




14.12 


14b Nd 




13.15 


natg m 


12.22 


13.86 


1 4S n 




15.01 






13.23 


1 ^4 i~t 

i54 Sm 




12.08 


™ a *Gd 


11.82 


13.06 


15b Gd 




14.15 


it>s Gd 




13.06 


i fin /^i i 
lbU Gd 




12.03 


" at Dy 


12.29 


12.97 






13.45 


1 f!4r> 

i04 Dy 




12.54 




13.04 


13.87 


lbb Er 




14.46 


1 

lb8 Er 




13.51 


i7U Er 




13.22 


Hi 


13.03 


15.13 


178xjr 

III 




1 C A A 

15.44 


Hi 




1 a on 

14.89 




12.36 


15.81 


182 W 




16.37 


184 W 




15.79 


186 W 




15.32 




12.74 


15.88 


198 Hg 




17.17 


200 Hg 




16.29 


202 Hg 




15.43 


204 Hg 




14.58 


nai Pb 


13.45 


15.15 


206p b 




15.54 


208p b 




14.97 


232 Th 


12.56 


13.71 


234 U 


13.79 


14.89 


236TJ 


13.09 2 


14.17 


238 |j 


12.57 2 


13.51 


242p u 


12.90 


13.13 


244p u 


12.40 3 


12.70 



TABLE I: Calculated rates for natural elements and for the important individual isotopes. All rates are in the 
form xx x 10 b s _1 . If the natural abundance has more than 90% of a given isotope then the Nuc column gives 
this particular isotope and the calculation is done for this nucleus only. Whenever a nucleus has the superscript 
nat the calculation is a weighted combination of all isotopes contributing more than 10% to the natural abundance. 
An empty experiment box means that the given isotope has not been measured. Natural abundance measurements 
appear in the box next to the most abundant isotope. Experimental data are from Ref. |4j, where the original 
sources may be found. 1 is taken from Ref. [24| , 2 is from Ref. [2l| and 3 is from Ref. [2{| 



